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Dynamics of a cold trapped ion in a Bose-Einstein condensate 
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We investigate the interaction of a laser-cooled trapped ion (Ba+ or Rb+) with an optically confined Rb 
Bose-Einstein condensate (EEC). The system features interesting dynamics of the ion and the atom cloud as 
determined by their collisions and their motion in their respective traps. Elastic as well as inelastic processes 
are observed and their respective cross sections are determined. We demonstrate that a single ion can be used to 
probe the density profile of an ultracold atom cloud. 
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The realization of a charged quantum gas, formed by laser- 
cooled, trapped ions and ultracold neutral atoms, offers in- 
triguing perspectives for a variety of novel experiments. In ad- 
dition to studying atom-ion collisions and chemical reactions 
in a regime where only one or few partial waves contribute, in- 
teresting phenomena such as charge transport in the ultracold 
domain llj], polaron-type physics |2-4], and novel atom-ion 
bound states ||5|] can be investigated. Further, the production 
of cold, charged molecules in a well-defined quantum state 
is an important goal in molecular physics [6, 7]. Techniques 
and tools that are widely used to control neutral atomic quan- 
tum gases, e.g. Feshbach resonances |8], may be employed 
in atom-ion collisions as well [9]. First observations of cold 
collisions with trapped atoms and ions in the mK regime have 
been made using Yb+ ions and a magneto-optical trap for Yb 
[10, 11] and very recently with a Bose-Einstein condensate 
(Yb^, Rb) lll2ll . Here we present our investigations where we 
use a defined small number of Ba+ or Rb+ ions which are 
immersed in an ultracold cloud of Rb atoms, which is either 
Bose-Einstein condensed or at sub-|jK temperature. We find 
that at our current settings micromotion of the trapped ions 
plays an important role in our experiments. The elastic cross 
sections are large as expected. Inelastic processes are in gen- 
eral strongly suppressed, which is important for planned fu- 
ture experiments. Charge transfer is directly identified in our 
experiments and found to be the dominant inelastic collision 
channel. As a first application we show how a trapped ion can 
be used to locally probe the density profile of a condensate. 

In the inhomogeneous electrical field of an ion, a neutral 
atom is polarized and attracted towards the ion. This inter- 
action can be expressed as a long-range polarization potential 
V{r) = —Ci^/lr'^ with C4 = q^a/Aneo, where q is the charge 
of the ion and a the dc polarizability of the atom. The char- 
acteristic radius of the potential is given by r* — ^/JiC^Jh, 
where ji is the reduced mass. For a ^^Rb atom (a = 4.7 x 
10^^' m^ 11411 ') interacting with a '^^Ba+ ion, this characteris- 
tic radius is r* = 295 nm. This is much larger than the typical 
length scale of the neutral atom-atom interaction potential as 
given by the van der Waals radius, which for ^^Rb is R^dw ~ 
4nm 1 8]. For collision energies £ above ^b x 100|jK, the cross 
section for elastic scattering can be approximated by the semi- 
classical estimate C7ei = 7t{nC4/Ji^y/^{l + 7t^/l6)E-^^^ iflll . 
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FIG. 1. Combined atom-ion trap: A linear radiofrequency Paul trap 
is used to store ions [20], whereas the ultracold atoms are confined in 
a crossed optical dipole trap. By precisely overlapping the positions 
of these two traps, a single ion can be immersed into the center of the 
ultracold neutral atom cloud. 



At lower energies a full quantum mechanical treatment is nec- 
essary, which takes into account individual partial waves. In 
addition to elastic collisions, inelastic processes can also oc- 
cur, such as charge transfer RbH-Ba+ — >Rb+H-Ba or molecule 
formation RbH-RbH-Ba+ ^RbH-(BaRb)+. 

Here we present the first results obtained with our novel hy- 
brid apparatus, where both species (Ba+, Rb) are first trapped 
and cooled in separate parts of a UHV chamber in order to 
minimize mutual disturbance and are then brought together 
We use a linear Paul trap into which we can load a well- 
defined number of ^^^Ba+ ions (typically 1) by photoioniz- 
ing neutral Ba atoms from a getter source. Photoionization 
is done with a diode laser at 413nm, which drives a reso- 
nant two-photon transition from the ground state to the con- 
tinuum via the ^Di state. The '^^Ba+ ion which has no hy- 
perfine structure is Doppler-cooled to mK temperatures us- 
ing a laser at 493 nm. Another laser at 650nm repumps from 



a metastable D3 m state. The ion's fluorescence is collected 
with a lens of NA=0.2 and then detected using an electron 
multiplying charged-coupled device (EMCDD) camera. The 



linear Paul trap (see Fig. 1) consists of four blade electrodes, 
which are placed at a distance of 2 mm from the trap axis, 
and two endcap electrodes at a distance of 7 mm from the 
trap center. It is operated with a radiofrequency (rf) drive 
of i2 = 27i:x5.24MHz with an amplitude of 1400Vpp and 
a dc endcap voltage of about lOOV. With these parameters 
we measure the axial and the radial trap frequencies to be 
coi,ax ~ 27rx80kHz and tOi,iad ~ 27rx200kHz, respectively. 

The sample of ultracold Rb atoms is prepared similarly 
as described in 111611 . Atoms from a magneto-optical trap 
are magnetically transferred over a distance of 43 cm into 
a quadrupole loffe configuration (QUIC) trap. Evapora- 
tive cooling yields a cloud of 3x10^ ^^Rb atoms in state 
I F = 1 , niF = — 1 ) at a temperature of about 1 |jK. From 
the magnetic trap we load the atom cloud into a vertical 1- 
dimensional optical lattice trap formed by two counterprop- 
agating laser beams (diameter w 500 |im) at a wavelength of 
1064nm with a total power of 2 W. We use the lattice as an 
elevator to transport the ultracold atoms upwards into the cen- 
ter of the Paul trap which is located 30 cm above the QUIC 
trap. The lattice is moved by changing the relative detuning 
of the two beams in a controlled manner II17I1 . The transport 
efficiency reaches more than 60%. During transport the depth 
of the lattice is sufficiently small (~ ^b x lOpK) so that evap- 
oration keeps the atom temperature at about 1 |jK. After trans- 
port the atomic sample is loaded into a crossed optical dipole 
trap (A = 1064nm). By lowering the trap depth, we evapo- 
ratively cool the atoms and end up with a BEC of about 10^ 
atoms confined in a dipole trap with trap frequencies of about 
(0^{, = (60 Hz, 60 Hz, 8Hz). The condensate is detected via 
standard absorption imaging. 

The BEC is initially located about 300 |jm away from the 
ion. We move the ion within 2 ms along the Paul trap 
axis into the BEC by changing the endcap voltages. The 
cooling lasers for the ion are switched off in order to en- 
sure that the ion relaxes into its electronic ground state 
\F = 1/2, mp = ±1/2) and to avoid changes in the atom-ion 
collision dynamics due to the cooling radiation. 

Initially the temperature of the ion is on the order of a few 
mK, which is much larger than the optical trap depth of about 
ksxl [iK. Thus almost any collision with a cold atom leads to 
the atom being lost from the trap. 

As a consequence, one could naively expect, that the atom 
loss stops after a few collisions, once the ion is sympatheti- 
cally cooled to atomic temperatures. However, when we mea- 
sure the number of Rb atoms remaining in the trap as a func- 
tion of the interaction time, we find a continuous loss of atoms 
(Fig . 2). It is the driven micromotion of the rf trap, which 
OlSIl is responsible for this continuing loss. In an atom-ion 
collision, energy can be redistributed among all motional de- 
grees of freedom, enabling also the flux of micromotion en- 
ergy to secular motion. After each collision, micromotion is 
quickly restored by the driving rf field. An equilibrium be- 
tween the energy that is inserted by the driving field and the 
energy taken away by the lost atoms is reached. Thus the min- 
imal temperature of a sympathetically-cooled ion stored in a 
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FIG. 2. Number of remaining Rb atoms as a function of time as 
a sample of Rb atoms interacts with an ion. (a) A single Ba+ ion 
is immersed into the center of a thermal Rb cloud. The line is an 
exponential decay fit. (b) A single ion (Ba+ or Rb+) is immersed 
in a Rb Bose-Einstein condensate. The lines are fits based on our 
simple model described in the text. 



rf trap is determined by the amount of micromotion of the ion 
|19J. Stray electric fields which shift the ion position away 
from the zero of the rf potential increase micromotion, lead- 
ing to the so-called excess micromotion BlSll . By applying dc 
electric fields along the radial direction this part of the micro- 
motion is tunable. 

Fig. 2 shows the elastic collision measurements with ei- 
ther (a) a thermal cloud of atoms (temperature rRb = 80nK, 
which is just above Tc) or (b) a BEC. In the following we 
analyze these data in detail starting out with the case for ther- 
mal atoms. In general, we have checked that losses due to 
collisions with the background gas are negligible. The loss 
of atoms is then described by A^ = —na^ivi, where n is the 
density of the atom cloud at the position of the ion, Cei the 
elastic atom-ion cross section and vi the velocity of the ion. In 
principal all three quantities are unknown. Based on our mea- 
surements and additional constraints of the following model 
we can still give estimates for them. 

It is important to note that the density profile of the atom 
cloud is depleted locally due to the collisions with the well- 



localized ion. According to its mass mj, energy Ei — miv^/2, 
and trap frequency coi the ion will be locaUzed inside a region 
(which for simplicity we choose to be spherical) of radius tq = 

\ Ei/mi(o\. Inside this sphere we assume a homogeneous 
density h. Right outside the sphere, the density is given by 
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FIG. 3. Loss rate of a thermal atom cloud when a static electric field 
is applied. In the presence of a field, the energy of the Ba+ ion and 
thus the amplitude of its secular motion are increased. As a result 
the sphere of depletion becomes larger, leading to an enhanced atom 
loss rate. For fields | (fdc l> 30 V/m, the amplitude of the secular 
motion is so large, that the ion spends a significant amount of time in 
a region of lower atom density. Therefore the loss rate decreases at 
high fields. 



n = Na)^^{mYib/'2KkBTjib)^'^, which is the peak density of a 
thermal cloud of Rb atoms (mass mRj,) confined in a harmonic 
trap with mean trap frequency dJRb = {co^^co^^co^^Y'^. The 
measured atom loss A^, which at f = is A^ = 1.5 x 10^ s^' 
(as can be read off from Fig. 2), must match the net flux of 
atoms into the "sphere of depletion", leading to the balance 
condition A^ = KrQ{n — «)vth. Here the thermal velocity of 
the atoms is given by Vth — ^ksT^h/^^ni^h- By setting n = 
in the balance condition we get a lower bound for the ion 
energy, which for our parameters is equal to A:bX 17mK. We 
can also estimate an upper limit of the ion energy from our 
method of compensating micromotion. The compensation is 
based on minimizing the position shift of the ion when the rf 
amplitude is changed. With this method we can reduce the 
DC electric field at the position of the ion to below 4 V/m, 
corresponding to a maximum ion energy of A:Bx40mK lIlSll . 
By taking the midpoint between the two bounds we estimate 
the ion energy to be about fc^xSOmK. Plugging this energy 
into the expressions above we get ro ~ 1 .45 pm and a density 
of n « 0.45 «. Moreover we can determine the cross section 
Cei to be 1.9 X 10^'"* m^. This value is in rough agreement 
with the semiclassical estimate cr^f™ = 9 x 10^'^m^. 



A slightly different analysis can be done for the measure- 
ment with the BEC shown in Fig. 2(b). Again, collisions of 
the atoms with the trapped ion will lead to a sphere of de- 
pletion within the condensate. However, the flux of atoms 
into the sphere is now driven by mean field pressure rather 
than thermal motion. From numerical calculations using a 
3D Gross-Pitaevskii equation with an absorptive (imaginary) 
potential term, we find that the flux into the sphere can be 
approximated by ^Ttr^hv^, as long as n > 0.1«. Here v^ = 
ATi^j2Ka{n — n)/m^\, is the velocity with which the atoms en- 
ter the sphere, a = 5.61 nm is the Rb-Rb s-wave scattering 



FIG. 4. Left: Fluorescence image of two Ba+ ions. Right: Fluores- 
cence of one Ba+ next to an unknown dark ion . We infer the 
existence of the dark ion from the position of the Ba+. 
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length and n = {l5m\^(ol^N /fia^l^) /2>n is the peak den- 
sity of the condensate. By equating the net flux to the mea- 
sured lossrate of the BEC at f = 0, A^ = 8 x 10"* s^^ (see Fig. 2), 
we find that the ion energy has to be larger than ^b x 1 mK. Our 
analysis shows, that on the other hand the ion energy has to be 
smaller than in the thermal case. Indeed, to minimize the tem- 
perature of the trapped ion, a comparatively large effort was 
made to compensate micromotion for the experiments with 
the BEC. Given that Oe\/olf™ « 2 as in the thermal case, our 
data suggests a reasonable ion energy of about A:b x 5 mK. This 
implies ro ~ 0.8 pm, a density of n « 0.1 n and a cross section 
of c^ei ~ 3.1 X 10^^'* m^. We also measured the loss rate for 
a single Rb+ ion and found it to be about the same as with 
Ba+ (Fig. 2). This is not surprising, since the atom loss rate 
depends only weakly on the ionic mass and the inner structure 
of the ion is irrelevant in the semiclassical regime. 

Following this discussion, we expect the net flux into the 
sphere and thus the atom loss rate to rise when we intention- 
ally increase the dc electric field. For small fields this is also 
what we observe in the experiment (Fig. 3). The maximum 
loss rate is observed at an electric field of (ode ~ 30 V/m, where 
the amplitude of the ion's secular motion is ro « 20|jm. This 
value is comparable to the size of the atom cloud, which has 
an extension of about 15 |jm along the radial and 80 pm along 
the axial direction. The model of a well-localized ion and a 
sphere of depletion is clearly no longer valid in this regime. 
We explain the decrease of the atom loss rate for even higher 
fields by the fact that the ion spends a significant amount of 
time in regions of lower atom density. 

In addition to the elastic processes discussed so far, we 
have also investigated inelastic atom-ion collisions. For this 
we load two '^^Ba+ ions into the ion trap. Typically after a 
time corresponding to lO'^-lO^ elastic atom-ion collisions, the 
fluorescence of one of the Ba+ ions is lost (Fig. 4). Since 
the position of the remaining bright Ba+ ion does not change, 
we infer that the other Ba+ ion has been replaced by an un- 
known dark ion formed in a reaction. We can determine the 
mass of the dark ion by measuring the radial center-of-mass 
(COM) trap frequency. For this we modulate the amplitude 
of the rf voltage. When the modulation frequency is equal to 
the trap frequency the ion motion is resonantly excited and 
a drop of the Ba+ fluorescence signal is observed. We only 
observe two resonance frequencies at 245 kHz and 400 kHz. 
We have checked that Ba+ gives rise to the 245 kHz reso- 
nance. Using the scaling properties of the analytic expression 
for the COM mode |20], the 400kHz resonance then corre- 
sponds to the mass of ^^Rb+. We conclude that the dominant 
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FIG. 5. Number of Rb atoms remaining in the trap depending on tiie 
position of tlie Rb+ ion relative to tlie center of tlie atom cloud. The 
measurement is performed with (a) a thermal cloud, (b) a partially 
condensed cloud and (c) an almost pure Bose-Einstein condensate. 
The interaction time was (a) 1.5 s, (b) 1 s and (c) 0.5 s. The solid lines 
are fits, where the ion energy Ej and the atom temperature Tri, are 
used as free fit parameters. 



inelastic process in our system is the charge transfer process 
Rb + Ba+ — >Rb++Ba with a cross section (Tch.ex. ranging be- 
tween 10^'^m^ and 10^^*^ m^. Our charge transfer results are 
comparable to the ones observed for the heteronuclear case of 
(Rb, Yb+) 1 13]. Homonuclear charge transfer rates are orders 
of magnitudes higher 1 1 1]. In our case the charge transfer is 
predicted to be dominantly radiative, where a photon carries 
away most of the 1 e V of energy released in this reaction 112 ill . 
The Rb+ ion can be stored in the trap for hours, even with- 
out any type of cooling, which opens up perspectives to study 
energetically resonant (Rb, Rb+) collisions. 

We now show how a single ion may be used to locally probe 
the atomic density distribution. By controlling the endcap 
voltage a single Rb+ ion is moved stepwise across the atom 
cloud and at each position x the atom loss is measured for a 
given interaction time (Fig. 5). The measurement is performed 
with various Rb clouds with different condensate fractions. 
We can theoretically reproduce the data shown in Fig. 5 with 
our model, which demonstrates our quantitative understand- 
ing of the dynamics. For this we write the total atom loss of 
a partly condensed cloud as a sum of the individual losses 
from the BEC and the thermal cloud as discussed before. 
For numerical calculations we choose cJei — Sx lO^^^m^ and 
make use of the well-known density profiles for a thermal 
and a condensed atom cloud confined in a harmonic trap. 
The thermal component and BEC are in equilibrium and obey 
Nc/N = 1 — {Tr\)/Tc)^ where Tc is the critical temperature and 
Nc/N is the condensate fraction which changes with time. The 
temperature Tr\, of the atomic sample is constant, as deter- 
mined by the optical trap depth. It is used as a free fit pa- 



rameter in our model together with the ion energy Ej, which 
we keep fixed for all three measurements. Fig. 5 depicts the 
total number of remaining particles N{x) as a function of the 
ion's position x. The values for T^b obtained from the fit (a) 
50 nK, (b) 35 nK and (c) 25 nK are in nice agreement with the 
temperatures determined separately in time-of-flight measure- 
ments. Moreover the fit suggests £/ « A:^ x 14 mK, which is 
in the same range as the temperatures found in the Ba+ ex- 
periments. The ion probe features a spatial resolution on the 
|im scale and has advantages compared to absorption imaging 
which integrates over the line of sight. 

In conclusion we have immersed cold trapped ions in a sea 
of ultracold neutral atoms in a novel hybrid apparatus with 
high spatial and temporal control. We have investigated the 
dynamics of the atom-ion system and have extracted first elas- 
tic and inelastic collision properties. In our present setup the 
collision energies are determined by the ionic excess micro- 
motion, which we plan to minimize for future experiments in 
the ultracold regime. 
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